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Fluctuations of scalar fields produced at the stage of preheating after inflation are so large that 
they can break supersymmetry much stronger than inflation itself. These fluctuations may lead to 
symmetry restoration along flat directions of the effective potential even in the theories where the 
usual high temperature corrections are exponentially suppressed. Our results show that nonthermal 
phase transitions after preheating may play a crucial role in the generation of the primordial baryon 
asymmetry by the Affleck-Dine mechanism. In particular, the baryon asymmetry may be generated 
at the very early stage of the evolution of the Universe, at the preheating era, and not when the 
Hubble parameter becomes of order the gravitino mass. 
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In the low-energy minimal supersymmctric standard 
model (MSSM) there exist a large number of D-flat di- 
rections along which squark, slepton and Higgs fields may 
get expectation values. We collectively denote them by tp. 
In flat space at zero temperature exact supersymmetry 
guarantees that the effective potential V(tp) along these 
D-flat directions vanishes to all orders in perturbation 
theory (besides the possible presence of nonrenormaliz- 
able terms in the superpotential) . In the commonly stud- 
ied supergravity scenario, supersymmetry breaking may 
take place in isolated hidden sectors jlj and then gets 
transferred to the other sectors by gravity. The typical 
curvature of D-flat directions resulting from this mecha- 
nism is M 2 ~ 1^-, where F is the vacuum expectation 

value (VEV) of the D-term breaking supersymmetry in 
the hidden sector. In order to generate soft masses of 
order of Mw in the matter sector, F is to be of order 
of (MvyMp) and sfermion and Higgs masses along flat 
directions turn out to be in the TeV range. 

D-flat directions are extremely important in the 
Affleck-Dine (AD) scenario for baryogenesis [|[ if large 
expectation values along flat vacua are present during the 
early stages of the evolving Universe. When calculating 
high temperature corrections to the effective potential at 
large tp one can often neglect the renormalizable interac- 
tions of the field tp with other fields, because the fields 
directly interacting with ip become too heavy to be ex- 
cited by thermal effects. Therefore thermal effects do not 
push the field <p> towards tp — in the early universe if 
initially it was large enough. 

For this reason, in the original version of the AD sce- 
nario HH it was assumed that the initial value of the 
scalar field tp was different in different parts of the uni- 
verse, as in the chaotic inflation scenario. However, re- 
cently it was pointed out that initial position of the field 
tp could be fixed by corrections to the effective potential 
which are proportional to the Hubble parameter H . In- 
deed, in generic supergravity theories moduli masses are 



of order of the Hubble parameter H [Q , || . One of the 
reasons is that the effective potential provides a nonzero 
energy density V ~ |P| 2 which breaks supersymmetry. 
If the vacuum energy dominates, the Hubble parameter 
is given by H 2 = (V/3Mp) and therefore the curvature 
along the D-flat directions becomes M 2 = off 2 , where 
c may be either positive or negative. (Here we use the 
reduce Planck mass Mp ~2x 10 18 GeV.) A similar con- 
tribution appears after inflation as well. For a mechanism 
of stronger supersymmetry violation see || . 

The most relevant part of the entire AD potential in 
presence of nonrenormalizable superpotentials of the type 
SW — (A/nM n ~ 3 )tp n somewhat schematically can be 
represented as follows [Q: 



V{tp)=c x m 2 3/2 \p\ 2 -c 2 H 2 \tp\ 2 + 
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AXHip n e lne + h.c. 
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Here Cj are some positive constants which we suppose to 
be 0(1), to 3 /2 = O(l) TeV is the gravitino mass, M is 
some large mass scale such as the GUT or Planck mass, n 
is some integer, which may take such values as 4, 6, etc., 
and the A-term violates the baryon or the lepton number 
and has a definite CP-violating phase relative to tp. This 
leads to symmetry breaking with tp <~ (H M n ~ 3 / A) 1 /™ -2 . 
For dcfinitcncss we will take here n — 4 and M = Mp. 

The main idea of ref. jj] is that in the early universe 
the field tp will be fixed at tp ~ y/HM P /X. When H 
drops down to m 3 / 2 , the terms proportional to H be- 
come smaller than the first term in ([I]), and the field tp 
begins oscillating near tp = 0. If the field tp can be asso- 
ciated with flat directions for squark-slepton fields, then 
the oscillations of the field tp near tp = in the presence 
of the A-term of eq. (^) may produce baryon asymmetry 
of the universe. 

The aim of this Letter is to show that there is a much 
stronger mechanism of supersymmetry breaking in the 
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early universe. This mechanism may lead to corrections 
to the quadratic terms in the effective potential of the AD 
field which are much greater than H 2 ip 2 and to A-terms 
much greater than the one of eq. ([!]). As a result, under 
certain conditions the Affleck-Dime mechanism of baryon 
asymmetry generation may become quite different from 
the version proposed in Q. 

The mechanism of supersymmetry breaking which we 
are going to discuss is related to particle production at 
the first stage of reheating after inflation. Indeed, Kof- 
man, Linde and Starobinsky have recently pointed out 
that a broad parametric resonance, which under certain 
conditions occurs soon after the end of inflation, may lead 
to a very rapid decay of the inflaton field [[tJ . The inflaton 
energy is released in the form of inflaton decay products, 
whose occupation number is extremely large. They have 
energies much smaller than the temperature that would 
have been obtained by an instantaneous conversion of the 
inflaton energy density into radiation. Since it requires 
several scattering times for the low-energy decay prod- 
ucts to form a thermal distribution, it is rather reason- 
able to consider the period in which most of the energy 
density of the Universe was in the form of the nonthermal 
quanta produced by inflaton decay as a separate cosmo- 
logical era, dubbed as preheating to distinguish it from 
the subsequent stages of particle decay and thermaliza- 
tion which can be described by the techniques developed 
in H . Several aspects of the theory of explosive reheat- 
ing have been studied in the case of slow-roll inflation ||] 
and first-order inflation ficf l and it has been recently pro- 
posed that non-thermal production and decay of Grand 
Unified Theory bosons at the intermediate stage of pre- 
heating may generate the observed baryon asymmetry 

& 

One of the most important consequences of the stage 
of preheating is the possibility of nonthermal phase tran- 
sitions with symmetry restoration jl2| , [[t3| . These phase 
transitions appear due to extremely strong quantum cor- 
rections induced by particles produced at the stage of 
preheating. What is crucial for our considerations is that 
parametric resonance is a phenomenon peculiar of par- 
ticles obeying Bose-Einstein statistics. Parametric res- 
onant decay into fermions is very inefficient because of 
Pauli's exclusion principle. This means that during the 
preheating period the Universe is populated exclusively 
by a huge number of soft bosons, and the occupation 
numbers of bosons and fermions belonging to the su- 
permultiplet coupled to the inflaton superfield are com- 
pletely unbalanced. Supersymmetry is then strongly bro- 
ken during the preheating era and large loop corrections 
may arise since the usual cancellation between diagrams 
involving bosons and fermions within the same super- 
multiplet is no longer operative. Therefore all results 
obtained in fL^| , |l3f| apply to the modification of the 
effective potential along the flat directions. As we will 
see, the curvature V"(ip) during the preheating era of- 
ten is much larger than the effective mass H 2 that L»-flat 
directions acquire in the inflationary stage. 



Let us consider, e.g., chaotic inflation scenario |TJ] 
where the inflaton field couples to a complex x-field, 

M 2 

V=^f<f> 2 + 9 2 cj ) 2 \ X 2 \ . (2) 

We take the inflaton mass ~ 10 13 GeV in order for 
for the density perturbations generated during the infla- 
tionary era to be consistent with COBE data. Inflation 
occurs during the slow rolling of the scalar field </> from 
its very large value. Then it oscillates with an initial am- 
plitude <fio ~ Mp. Within a dozen oscillations the initial 
energy ~ 0q is transferred through the interaction 
g 2 <p 2 x 2 to bosonic %-quanta in the regime of parametric 
resonance JjJ. At the end of the broad parametric res- 
onance the field cf> drops down to <fi e ~ 10 17 GeV. An 
exact number depends logarithmically on coupling con- 
stants; we will use 4> e ~ 10 17 GeV ~5x 10 -2 Mp for our 
estimates. After this stage, the universe is expected to be 
filled up with x-bosons with very large occupation num- 
bers rik ~ .9 2 , with relatively small energy per particle, 
E x ~ y/ gM<p(f) e ~ Q.2yJ gM^M-p. It is especially impor- 
tant that the amplitude of field fluctuations produced at 
that stage is very large , Jl2| , 

(X 2 > ~ 5 x ICTV^Mp . (3) 

If, for example, the energy of the inflaton field after pre- 
heating were instantly thermalized, we would obtain a 
much smaller value (x 2 ) ~ 10~ 4 M < pMp. In realistic mod- 
els thermalization typically takes a lot of time, and the 
value of (x 2 ) after complete thermalization is many or- 
ders of magnitude smaller than 

Note, that large fluctuations (||) occur only in the 
bosonic sector of the theory, thus breaking supersymme- 
try at the quantum level. Anomalously large fluctuations 
of (x 2 ) lead to specific nonthermal phase transitions in 
the early universe |L2j , (l3| . It will be interesting to study 
possible consequences of this effect for supersymmetric 
theories with flat directions. 

Let us make the natural assumption that the field tp 
labelling the Z?-flat direction couples to the generic super- 
multiplet x by a renormalizable interaction of the form 
£int = ^ 2 M 2 |x| 2 coming from some _F-term in the po- 
tential. Reheating gives an additional contribution to the 
effective mass along the (^-direction : 

AM 2 ~ 2 h 2 (x 2 ) ~ KT 1 j M<p M P . (4) 

As a result, the curvature of the effective potential be- 
comes large and positive, and symmetry rapidly restores 
for 1CT 1 Mcf, M P > c 2 H 2 . At the end of preheating 
in our model H ~ 2 x 1(T 2 M . Thus, for c 2 = 0(1) 
symmetry becomes restored for h 2 > 2 x 10~ 8 g. In 
other words, supersymmetry breaking due to preheating 
is much greater than the supersymmetry breaking pro- 
portional to H unless the coupling constant h 2 is anoma- 
lously small. This leads to symmetry restoration in the 
AD model soon after preheating. 
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a baryon asymmetry ub/ti^ = O(l). 

In our case baryon number production occurs at very 
early stages, during the preheating era. This is possi- 
ble because the same nonthermal effects generating a 
large curvature at the origin give rise also to sizable 
CP-violating terms along the AD flat direction. In- 
deed, under very general assumptions, one can expect 
the presence in the Lagrangian of nonrenormalizable CP- 
violating terms of the type afy 2 /?^™ - )ip m +h..c., where 
m > 2. During the preheating period large amplitudes 
of the x-field give rise to the CP- violating term 



Let us check the applicability of our results. First of 
all, parametric resonance takes place for gcf> > M^, where 
<j> ~ 10 _2 Mp at the end of preheating. This implies that 
g > 1CP 4 . Secondly, even though self-interactions of the 
X-field do not terminate the resonance effect since parti- 
cles remain inside the resonance shell, creation of quanta 
different from x may remove the decay products of the 
inflaton away from the resonance shell. This leads us to 
consider two different possibilities. The resonance does 
not stop if scatterings are suppressed by kinematical rea- 
sons, i.e. if the non-thermal plasma mass of the final 
states is larger than the initial energy of the x's. If the 
final states are identified with the ^-quanta, this hap- 
pens for h > g. Otherwise, scatterings occur, but are 
slow enough to not terminate the resonance if the in- 
teraction rate T ~ n x er, where a denotes the scattering p or m _ 3 we g e ^ a CP-violating term 
cross section, is smaller than the typical frequency of os- 
cillations at the end of the preheating stage ~ gix 2 ) 1 ^ 2 - 
Identifying again the final states with the </?'s, this condi- 
tion translates into h 2 < 10 2 g, which is compatible with 
the condition for symmetry restoration discussed above. 

There is another necessary condition in order to have 
broad resonance and a rapid production of x-particles: 
the constant contribution h\tp\ to effective mass of the 
field x induced by the condensate \ipo\ should be smaller 



(X 2 



f m e imti 



h.c. 



(5) 



h.c. 



than the typical energy of the decay products E x . This 
translates into a bound on the coupling A when taking 
into account the constraint on h discussed previously. 
The exact bound depends upon the choice of the index 
n and the mass M. For n = 4 and M = Mp one has 
the condition h 2 < 3gX. This condition simultaneously 
guarantees that the energy density of the AD field ip is 
smaller than the inflaton energy density at the end of 
the broad resonance. All these constraints are not very 
restrictive. For example, one may take h = 2g to satisfy 
the condition h > g. Then one has the constraints g > 
10" 8 and A > g. 

Let us now describe the dynamics of the yj-field during 
preheating. Under the conditions described above, the 
AD field if rapidly rolls towards the origin ip = and 
makes fast oscillations about it with an initial amplitude 
given by f ~ \fo\- The frequency of the oscillations is of 
order of (AM 2 ) 1 / 2 and is much higher than the Hubble 
parameter at preheating ~ 10 _1 M$. The field is under- 
damped and relaxes to the origin after a few oscillations. 

The crucial observation we would like to make here is 
that baryon asymmetry can be efficiently produced dur- 
ing this fast relaxation of the field to the origin. Our 
picture is similar to that of ref. , but there are consider- 
able differences. According to B, the baryon asymmetry 
is produced at very late times when the Hubble param- 
eter H becomes of order of the gravitino mass TO3/2 ~ 
TeV. At this epoch the field becomes underdamped and 
the CP- violating A-terms in the Lagrangian are com- 
parable to the baryon-conserving ones. The field feels a 
torque from the A-terms and spirals from the initial point 
if = \fo\e 10 inward in the harmonic potential, producing 



with A ~ 1CT 1 g~ l a M^. 

Note that fluctuations (x 2 ) were absent during infla- 
tion; they increase very sharply (exponentially) during 
preheating, and therefore they begin strongly affecting 
the position of the scalar field only at the end of the 
broad parametric resonance, when they reach their max- 
imal value. Thus, to the first approximation one has the 
following initial conditions fo r the field f at the end of 
the broad resonance: \ip \ ~ \J HMp /A ~ ^/A/^Mp/lOA, 
ipo = 0. As for the complex phase 9, its initial value is 
determined by the P-dependent A-term ([!]). In general, 
this phase is quite different from the phase determined 
by eq. (5). As a result, the field tp will spiral down to 
ip = acquiring baryon charge density ns — 2\f 2 \9. 

After the end of the stage of broad parametric res- 
onance the value of (x 2 ) decreases approximately as 
a~ 2 ~ t~ l . As a result, both AM V and A after preheat- 
ing decrease as t~ x , i.e. in the same way as the Hubble 
constant H in the radiation dominated universe after pre- 
heating. Therefore the equation describing the relaxation 
of the ip-field toward the origin 



if + 3 H if ■ 



AM 2 ip - 



A if*) 



takes the following form: 











(6) 



(7) 



Before going further, let us compare the magnitudes of 
the last two terms for \ f 1 ~ y/HM P /A - ^M^Afp/lOA. 
One can easily see that the P-violating term is greater 
than the mass term for a > 10 3 /i 2 \/A. If, e.g., one takes 
h ~ 10~ 7 , A ~ 10~ 2 , one finds that baryon number vio- 
lation may be substantial even if the coefficient a in the 
A-term (||) is extremely small, a ~ 10~ 12 . 

It is convenient to introduce new variables y — 
fyf Xt/Mp and r = h^/2tMp/g. In these variables eq. 
(0) simplifies: 



y + y + 



aV2 r2 



t h yfgX 



(8) 
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In the new variables the motion of the field y = \y\e %s 
begins at \y$\ = 1, at the initial moment To ~ 10 4 h/y/g. 
The condition of vanishing initial velocity of the field ip 
translates into the condition y' = yo/ro. To calculate the 
ratio tib/s, where tib — 2\{p\ 2 Q 1 and s is the entropy den- 
sity, we will introduce a fictitious reheating temperature 

t r ~ io"V m * m p ~ 5 x 1013 GeV - This is the tem " 

perature which would be reached by our system if ther- 
malization would occur instantaneously after the end of 
the broad parametric resonance at 4> e ~ 10 _1 M P . Even 
though thcrmalization may occur much later, this con- 
cept may be quite useful because at the radiation domi- 
nated stage the energy density of the universe decreases 
in such a way that at the moment when thermalization 
actually occurs the resulting temperature T will be equal 
to the redshifted value of the fictitious reheating temper- 
ature, T ~ Tr^/tq/t. This leads to the following ex- 
pression for the baryon number B = ur/s soon after 
preheating: 



models. We have found that if the inflaton field (j> couples 
in a renormalizable way to x-bosons, which are weakly 
coupled to the AD field ip, then the effect of parametric 
resonance may induce very large masses for particles cor- 
responding to flat directions of the AD potential. The 
same effect may induce large terms violating baryon con- 
servation (|J). As a result, one can obtain large baryon 
asymmetry even in the models where all relevant coupling 
constants are extremely small. 

We would like to thank D. Lyth, LA. Kofman and E. 
Kolb for many useful discussions. GA is supported by 
the DOE under contract DE-AC02-76CH03000; AR is 
supported by the DOE and NASA under Grant NAG5- 
2788; AL is supported in part by the NSF grants PHY- 
9219345 and AST-9529225. 
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102 T R 



2 x 10~ 2 |2/| 2 f 



AV5 



(9) 



We solved eq. (0) numerically for various values of pa- 
rameters and calculated the ratio of n B to the entropy 
density s. We have found that this ratio oscillates and 



approaches a constant at large times. For 



< 1 and 



h > 10 _4 v / g the typical value of the baryon asymmetry 
is given by 



B = ^ 



1(T 



/(to) ■ 



(10) 



Here /(to) is a certain function of to ~ 10 4 /i/ v / 5 : /(l) ~ 
1, /(10) ~ 0.1, /(TOO) ~ 0.05. 

Validity of this result depends on details of the theory; 
typically it gives a reliable estimate of the baryon asym- 
metry only for B <C 1 (which is what we need), but even 
in this case some care should be taken. For example, one 
can show that for h 2 < gX the temperature T remains 
much greater than the time-dependent effective mass of 
the field ip until this mass approaches the constant value 
0(7713/2). In this regime quarks acquire large effective 
mass 0(T), and the field ip cannot decay and transfer its 
baryon asymmetry to fermions until temperature drops 
to 0(7713/2). This may lead to some corrections to eq. 
(0) Hi • Also, the simple scaling rules used in the deriva- 
tion of eq. (|1C| ) may break if, e.g., at some intermediate 
stage the universe becomes dominated by nonrelativistic 
particles. 

A detailed investigation of baryon asymmetry produc- 
tion in realistic theories including all of the effects men- 
tioned above should become a subject of a separate in- 
vestigation. The main purpose of our paper was to show 
that parametric resonance and nonthermal phase transi- 
tions found in J? 12 may lead to strong supersymmetry 
breaking in the early universe and to considerable modi- 
fications in the theory of baryogenesis in supersymmetric 
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